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We study the charge and spin currents passing through a magnetic tunnel junction (MTJ) on 
the basis of a tight-binding model. The currents are evaluated perturbatively with respect to 
the tunnel Hamiltonian. The charge current has the form A[Mi(t) x Mi(t)] ■ Mj + BMi{t) ■ M2, 
where Mi{t) and M2 denote the directions of the magnetization in the free layer and fixed 
layer, respectively. The constant A vanishes when one or both layers are insulators, while the 
constant B disappears when both layers are insulators or the same ferromagnets. The first term 
in the expression for charge current represents dissipation driven by the effective electric field 
induced by the dynamic magnetization. In addition, from an investigation of the spin current, 
we obtain the microscopic expression for the enhanced Gilbert damping constant Aa. We 
show that Aa is proportional to the tunnel conductance and depends on the bias voltage. 
KEYWORDS: spintronics, magnetic tunnel junction, spin current, spin dynamics 



1. Introduction 

Magnetic tunnel junctions (MTJs), which consist of a thin tunnel barrier sandwiched be- 
tween two ferromagnetic layers,'"^ are promising for their use in magnetic random access 
memory (MRAM).^ However, the primary disadvantage of conventional MRAM designs, 
which employ a current-induced field to write data, is that the writing current increases 
with the device density. Thus, there has been considerable interest in exploiting spin-transfer 
torque (STT)^'^ instead. ^"^"^ In such an STT MRAM device, the critical current is proportional 
to the product of the volume and the Gilbert damping constant a of the free layer, making 
low a an important criterion for electrode materials. 

To this end, several studies have explored the dynamics and the distribution of the magne- 
tizations in STT MRAM by using the Landau-Lifshitz-Gilbert (LLG) equation with an STT 
term.^"^ However, other torques (spin torques) also act on the dynamic magnetization in the 
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free layer, which form in reaction to the outward flow of spins from the layer: Mizukami et al. 
experimentally showed that a increases with the thickness of the nonmagnetic metal (NM) 
layer in NM/Py/NM films, and that this enhancement continues up to thicknesses of several 
hundred nanometers. ^'^ Their experiment supports the importance of spin torques in the mag- 
netization dynamics of mesoscopic devices such as STT MRAMs. Further, this experimental 
finding was supported immediately by Tserkovnyak et al. 's^°'^' theory of spin pumping based 
on scattering theory, with additional theoretical confirmation by Umetsu et al. on the basis of 
the Kubo formula.^^'^^ 

Several studies have also investigated charge transport in the presence of magnetization 
dynamics in magnetic multilayers. It is known that dynamic magnetizations induce an efl'ec- 
tive electromagnetic field.^"*'^^ Ohe et al. simulated the effective electric field induced by the 
motion of the magnetic vortex core in a magnetic disk,^^ and the field was observed experi- 
mentally.^' Furthermore, Zhang et al. phenomenologically derived the LLG equation having 
the STT term induced by this eff'ective electric field. And Moriyama et al. observed the dc 
voltage across generated by the precession of the magnetization in an Al/A10v/Ni8oFe2o/Cu 
tunnel junction. The origins of this voltage have been discussed from a theoretical stand- 
point (scattering theory). ^''"^^ In addition, charge and spin currents in ferromagnets with mag- 
netizations that slowly vary in space and time have been studied microscopically.^'' ''^ These 
studies employed the s-d model in continuous space and treated the perturbation within the 
framework of the Keldysh-Green function.''^'''' 

Similarly, our aim is to describe the charge and spin transport in MTJs in the presence of a 
voltage across the barrier and the dynamical magnetization in the free layer. This situation just 
corresponds to an STT MRAM cell during the writing stage. In this paper, we microscopically 
describe the charge and spin currents passing through an MTJ. However, in contrast with 
previous works that relied on models in continuous space, we calculate the currents on the 
basis of a tight-binding scheme. This makes it easier to account for the properties of materials 
and the space dependence of the magnetization in magnetic multilayers, such as MTJs, with 
strongly inhomogeneous magnetic structures. In the calculations, we consider the voltage and 
the dynamics of the magnetization in Berry's adiabatic approximation under the assumption 
that the eff'ective exchange field is larger than the voltage and dynamics. Our model shows that 
the charge current induced by the dynamical magnetization has the form A[Ml(0 x Mi^{t)] ■ 
Mr + BMiSt) ■ Mr, where Ml(0 and Mr denote the directions of the magnetization in the 
free layer and fixed layer, respectively. The first term tends to the form given by Tserkovnyak 
et al.,^^ which expressed the dc current due to the precession of Ml(0 about Mr as a special 
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Fig. 1. Schematic of one-dimensional magnetic tunnel junction. Tlr is the tunneling amplitude and tij repre- 
sents the hopping matrix between sites / and j located at either side of the interface. MlCO and Mr denote the 
directions of the effective exchange fields for the left (L) and right (R) hand side layer, respectively. 



case; in this sense, our result is a generalization of their work. Furthermore, from the results 
concerning spin transport, we successfully derive the enhanced Gilbert damping and propose 
a microscopic expression for it. 

2. Model and Formalism 

2. 1 Model Hamiltonian 

We consider the motion of electrons in an effective exchange field. Furthermore, assume 
that the ferromagnetic layer on the left-hand side (LHS) of the MTJ is the free layer; that 
is, the direction of the field at time t in this layer, Mi^{t), rotates time-dependently (see Fig. 
1). Thus, the direction of the field on the right-hand side (RHS) (fixed layer), Mr, is time- 
independent. Note that we ignore the inner structure of the tunnel barrier and account for 
its properties via the simple tunnel amplitude Tlr between sites L and R, which denote the 
surfaces on the LHS and RHS, respectively. In this model, the total Hamiltonian for the MTJ 
is the sum of the one dimensional tight-binding Hamiltonians in the ferromagnetic layers, 

'Hl(0:= Yj A[-^ij'^-5ijJi.Mdt)-d]cj, (1) 

ijeLHS 

'Hr:= Yj c][-tifl-6ijJ^M^-&\cj, (2) 
and the tunnel Hamiltonian, 

•Ht := -Tlrc^cr + H.c, (3) 



J. Phys. Soc. Jpn. 



FULL PAPERS 



where cJ^(c,o-) is an operator that creates (annihilates) the a spin electron at site /, and tij is 
the hopping integral between sites i and j. The constant Ji^(Jr) represents the strength of the 
interaction between the spin of an electron and the effective exchange field on the LHS (RHS) 
layer; and & is the Pauli matrix, where hat denotes a 2 x 2 matrix in spin-space. 

2.2 Adiabatic approximation 

Assuming yL ^ ^|dML(0/d4 we adopt Berry's adiabatic approximation^^ for 'Hi^{t): 

Ciit) ^ Ui^{t)e^^'^^" di for / e LHS, (4) 

■^^0^^':= Yu d][-Uj\-6ijJi^&^dj, (5) 

!,;eLHS 

where c,(0 is in the Heisenberg representation with respect to 'Hi^(t), Ui^{t) is a rotation matrix 
satisfying the equation Ul^(t)Mi^(t) ■ &Ui^(t) = and yit) is Berry's phase defined by 



7(0 := i Jd? 



uUt) 



(6) 

TT 



With the approximation (4), we replace "Kt with 

^fit) := -TLR^e-^^W'^ Ulm^d^ + H.c, (7) 

where Ur is a rotation matrix satisfying the equation ^/^Mr • o-Ur = (T", and J, := 
f/^c, for i e RHS. Finally, our total Hamiltonian is 'Hit) := + "Hr + 'Hfit), where 
'^R = Z/,yeRHs \-tijl - %7R(T^j dj. Thus, a nonequilibrium statistical average of the form 
(dia-(t)dl^,(0) can be derived perturbatively with respect to 'J-{!f(t) using the Keldysh-Green 
function technique. 

2.3 Charge and spin currents 

The charge current P{t) and spin current P{t) passing through the MTJ are defined by 

nt):=2%jTRi^{dlmlUdt)^^^'^'''dL(t)) [1/s], (8) 

Pit) := 2%^Tr^ {dlit)Ul&Ui^(t)^^^'^'^dL(t)) [1/s], (9) 

where (• • •) denotes a statistical average in 'Hif)?^'^'' 
By introducing the lesser function, 

[^lr(^0]_, := ^{[4(O^y^,[t/L(0e'^«*=JL(0]J, 
eqs. (8) and (9) can be written in the form 

I\t) = 2'KTR^ixGlR{t,t), (10) 
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(11) 



(12) 



(13) 



(14) 



(15) 



In the first order in "H^it), we have 

Gl^it, t) - - Jlr J df't/L(Oe'^^'^'^fL(? - t')t-''''^^''Ul{,t)A{t, t'W^g^it' - tWl 
where < denotes the lesser component of Keldysh-Green functions,^''' and 

A{t,t') := f/L(Oe'^^'^"'e-^^('')*'[/^(f'). 
Moreover, we introduce the unperturbed Keldysh-Green functions defined by 
[gdt)U> := -^{T:dUt)dl,\ = -'-{TdUt)dl\6^^., 

where T is the time-ordering operator on the Keldysh contour, and (■ • ■)o denotes an equilib- 
rium statistical average in 'H^ + "Hr. Since IlCO is the diagonal matrix in spin-space, gi^{t) 
and Berry's phase factors commute. Thus, G^^S^, t) reduces to 

2;^e-'^''^'^LWM^)f^L(0A(r,£')[/;^R(£-£')t/R , (16) 

where we employ the Fourier transform of a function f{t) with respect to t, defined by the 
relation 

f{E) := J dre'^'^VCO. (17) 

E' in eq. (16) represents the energy that an electron obtains from the dynamics of the magne- 
tization; we consider it in the first order: 



GlM,t)^-T, 



LR 



{-ul 



{t)gaE)ULmigj,{E)u^ 



-T^LR 

Then, using the relations 



r d£ . n dA(t,t') 



df 



t/«f/R 



dE 



n dA{t, f) 



i dt' 



t'=t 



n dML(r) 



Ulitrg^EWdt) = IME) + d- ■ MdOAgdE), 



UlMEW^ = IME) + d- ■ M^Ag^(E), 



where 



gL(R)(E) := -tr|L(R)(^) = 



[gUR)(E)U + [guR)(E)] 



(18) 

(19) 

(20) 
(21) 

(22) 
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1 , [tL(R)(£)]t. - [gUR)(.E)].. 

Agu^m := -tr(t^|L(R)(£) = 



we can decompose G^^{t, t) into two terms 



Gl^{t,t) = \Gl^it,t) + &-GtJ.t,t). 



(23) 



(24) 



Here we define 



GfJj^, t) :- -Tlr 



d£ 



2n% 



gL(.E)ME) + Mdt) ■ MKAg^{E)Ag^{E) 



+ -Ml(0 X —r— ■ MRgL(£) -7^ It ,^ ■ MrZI^lC^) TT^ 

2 d? d£ 2 at aE 



GUt,t) :=-ri 



LR 



dE 



Infi 



(25) 



Mdt)Agi^{.E)g^{E) + MRME)AgKiE) + iM^t) x MRAg^{E)Ag^{E) 



, ^^^^ dML(0 _ ,^, dME) .ndMdt). ,j,,dUE) 



+ 1- { Ml(0 X 



dt 



^ - rr.. ^SR(.E) ^ fl dM^t) ^ . , . ,j,,dAg^(.E) 

X Mr§l(£) — TP — + o — 7— ^ M^Agi^iE) 



dE 



2 dt 



dE 



F(t) and P(T) axe expressed in terms of G^^(t, t) and G^^{t, t) as follows: 

r{t) = A'KT^Gl^itj), 
I\t) = A%T^i^Gl^{t,t). 



(26) 

(27) 
(28) 



Finally, taking the lesser components, we obtain the following in the low-temperature limit: 



r{t) = 2n\Ti^RnpL(juWR(M)Mdt) x 



dE 



dt 



Mr 



Api^iE) — — — Apji(E) 



dMdt) 



-r 

m = \ J dE[ApdE)AxR(E) + AxdE)Ap^(E)]Mdt) x Mr 



+ 



4nm 



LRl 



Pi^(jj)pR(M)-Mdt) X 



r 



dE 



dE 



A dAxL(.E) _ 
Ap^iE)—- — — pr(£) 



dE 



dE 



2 dt 



- rr.. ^XR(.E) dxdE) 

Pl(E) — -—Apr(E) 



dE 



dE 



Mdt) X 



dML(0 
dt 



X Mr 



+ Ap^(ju)ApM^^^^ X Mr 
2 d? 



(29) 



(30) 



^ 11 ^ 
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Here // is the chemical potential of the system. pl(r)(£') and zJpl(R)(£') are the spin- averaged 
local density of states (DOS) and the spin polarization of the local DOS, respectively, at the 
LHS (RHS) layer surface, defined by 

Pl(r)(£) := ~3glR)(E), (31) 

ApuR)iE) := -^3zjgL(R)(£), (32) 

where g^'s are retarded Green's functions from the calculations taking the lesser component. 
Furthermore, the;^f's are defined as the real parts of the retarded Green's functions, 

XuR)(E) := -^^L(R)(^). (33) 
n 

AxuKm ■= -^^8uR)iE). (34) 
n 

3. Discussion and Summary 

3.1 Charge current 

The form A[Ml(0 x Ml(0] • Mr + 5Ml(0 • Mr of the charge current (eq. (29)) driven 
by the magnetization dynamics is consistent with previous works; The first term tends to the 
form given by Tserkovnyak et al?^ in the special case where Ml(0 precesses around Mr, 
as discussed in §3.3. And Xiao et al.^'^ have derived the same form for the charge current 
passing through the MTJ on the basis of scattering theory in the continuum space, whereas, 
we calculate the current on the basis of the tight-binding model. New insights which we found 
out in this study are as follows. If the electronic structure in the two layers of the MTJ is the 
same [i.e., Api^{E) = Ap^{E) and Axl(E) = A^RiE)] or both layers are insulators, we have 
B = 0. However, if either layer is metallic, finite B should be measured because the real 
part of the retarded Green's function remains finite, which reflects virtual transitions through 
forbidden bands. The term [Ml(OxMl(0] "Mr in eq. (29) represents the charge current driven 
by the effective electric field (i.e., the spin electric field), as mentioned in §3.3. In other words, 
the effective electrochemical potential of the free layer is changed by the dynamics of Ml(0> 
and the resultant difference in electrochemical potentials between the two layers manifests as 
a bias voltage.^' This situation may be realized when a barrier exists between the electrode 
and lead, or when the diffusion constant of the free layer is small enough to maintain the 
changed effective chemical potential. Otherwise, this charge current will flow back to the 
reservoir connected to the free layer without tunneling through the barrier of the MTJ. 
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3.2 Spin current 

The term Ml(0 x Mr in eq. (30) represents the static effective Heisenberg coupling be- 
tween Ml(0 and Mr. That is, the equation of motion for Ml(0 described by this spin current 
corresponds to the equation = ^^^M^it) x Mr [5l(0 is defined by eq. (42)]. This 

aff"ords a Heisenberg coupling energy of -J^ffMi^it) ■ Mr, where 

7eff := -4;r|rLRp J' dE [Ap^iE)AxK(E) + Axl(E)Ap^(E)] (35) 
" dE3GlUE)A^(E)G'^, (E)AdE), (36) 



if 



G\^{E):=g]^{E)Tijg)^{E\ (37) 

ME) ■= g]^{Er' - g]^{Er\ (38) 

A,^^) describes the exchange splitting at site /, and this result agrees with the expression 

presented by Liechtenstein et al?'^ 

Let us consider the term Ml x in eq. (30): 

, dML(0 n - dML(0 

2n\T^^?pd^DpMMl^{t) X — ^ = -^rML(0 X — (39) 

at 2e^ at 

where e > is the elementary charge and f is the tunnel conductance of the MTJ, 

r := pMpr(m). (40) 

n 

This term describes the spin pumping in the MTJ and affords the following microscopic 
expression for the enhanced Gilbert damping constant: 

n r 

Aa = — , (41) 

where 5l(0 is the total spin polarization of the electrons in the LHS layer, 

Sl(0:=2^ r dEApiiE)M^it). (42) 

(SLHS 

Equation (41) agrees with the corrected Gilbert damping constant derived by Zhang et al?^ 
phenomenologically after considering the effect of the spin electric field induced by the dy- 
namic magnetization. In addition, in the present formulation, from the fact that Aa vanishes 
if one ignores Berry's phase (6),"^" it follows that one of the origins of spin pumping is the 
spin electric field. As a consequence of this, Aa is proportional to the conductance f . 
The size dependence of Aa can be described as follows: 

Aa oc\, (43) 
A 

where A is thickness of the free layer, because |5l(0I is roughly proportional to the volume of 
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the free layer, and T to the cross-sectional area of the barrier. 



3.3 Analysis of effective field 

For a more transparent physical interpretation of the currents, we rewrite eqs. (29) and 
(30) as follows: 



-eP{t) 



o-=±l 



(44) 



-em = 



eJeff fi dMr (0 



xMr 



+ 



n ^ 2e dt 



o-=±l 



where the "conductances" are defined by 
2;r|rLRpg^PLO")PRtrO^) 

n 

47r\T^^\^e^JpM^pR(M) 



AY : = 



7o 



n 



2n\T^^\^e' 



n 



2n[r^\'e' 

n 



dE 



dE 



Put(E) — — — — Apr(E) 



dE 



dE 



PrAE) — -= -= — Api^iE) 



dE 



dE 



and the effective driving fields can be defined by 



1 . . crfi _ dMdt) 



4(0 : = 



crfi dML(t) 



(45) 

(46) 
(47) 
(48) 

(49) 

(50) 
(51) 



2e dt 

The conductances represented by a capital letter denote the "Fermi surface terms," whereas 
those represented by a small letter denote the "Fermi sea terms." The spin-dependent effective 
voltage e^(0 • Mr in eq. (44) just corresponds to the spin electric field between the layers. 
To compare the expressions obtained in continuous space and in discrete space, let us define 
the correspondences M(r, t) := Ml(0 and M(r -I- Ar, t) := Mr, where Ar denotes the barrier 
thickness. Then we find s^t) ■ Mr ^ Ar' (-^) x ■ M(r, t), which is well-known as 
the spin electric field. When Ml(0 steadily precess about the direction of Mr with a constant 
cone angle 6 and a constant frequency a», the voltage is time-independent: 



4(0 ■ Mr = -a— sm' e, 

2e 



(52) 



This affords an estimate ftoj/le ~ 20 yuV at 10 GHz. The Fermi sea term in eq. (44) vanishes 
in this case. This result is in good agreement with that of Xiao et al.^'^ and Tserkovnyak et 
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al.^^ Note that in general the Fermi sea term is certainly the ac current. 

Next, let us consider the spin current (45). The terms including r^e^(0 + |^y);e^(0 • 
^^r]^l(0 describe the spin transport due to the spin cr component of the charge current. 
By considering as a driving force, we can interpret the term + Ml(0 • MRyJ^j e^(0 
as the "tunneling magnetoresistance (TMR) effect" in spin transport. 

3.4 Effects of bias voltage 

Finally, we consider the charge and spin transport in the presence of a bias voltage V{t) 
across the MTJ. In Berry's adiabatic approximation under the assumption yL ^ e|V^(OL the 
effects of V(t) can be included by replacing eq. (4) with 

c,a)^e-^/*^'«f/j^(^)eirW'^j. for / e LHS. (53) 
In the first order in the effective exchange constant and conductances differ as follows: 



Jes ^ ^eff + (7t - yi)-y(t) + ^r^^ In 



dt ' 



(54) 



AY ^AYll In 

dyu 



PRcriu) 

^Pl(m) 

^pM 



eV{t) 



eV{t)- J' 



aE 



d£3 



dE 



pRaiE) endV(t)] 
Pl(m)Prc-(m) 2 dt j 

Ap^E)'^^ - t^Ap^iE) ehdVit) 



, (55) 



Api^(ju)ApR(p) 



2 dt 



(56) 



7o 



n 



dE 



PlAE) 



d^AxRiE) , d^XLAE) 



d£2 



d£2 



^PRiE) 



eV{t) 



+ 



27:\T^R?e^ 



n 



dpuri/j.) dAp^ip.) dVcrO") . d'^ApRip.) 
2 — ^Pr(/^) -Pl<t(m) 



dfi dfi 



d//2 



dyu2 



eh dV{t) 
y dt ' 



(57) 



n 



r 



dE 



,^, d'AxL{E) , d^XRAE) , 
PrAE) — — + — TT. — ^PhiE) 



d£2 



d£2 



eV(t) 



+ 



2n\T^R?e^ 



n 



dpR^iju) dAp^Qu) d^pR^ijj) d^ApLip.) 

^ APl(m) -pRcriM) 



eh dVit) 



dfi dfj. dfi^ r^^' ,..^v-. 

In addition, a term describing the TMR effect, 

1 r- 1 1 r H fidV{t) 

- r + AYMAt) ■ Mr V{t) + — r + AyMi^it) ■ Mr -— ^ 
e ^ ^ -e 2 dt 

appears in the charge current, where 



2 dt 



(58) 



(59) 



n 



- ,^, d2;pR(£) , d^xAE) _ 
pAE) + -—^pAE) 



d£2 



d£2 



4^|V 

n 



Ay : = 



dE 



A .^. d2j;rR(£) , d^xAE) , 
ApAE) — —i — + — — 5 — ApAE) 



d£2 



d£2 
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For the spin current, a term describing the STT effect, 
1 



e 



(rL-rL)y(o-(r^ + 7f) 



MlCO + - 
e 



2e d? 



is added, where 



27r|rLRpeVLtrOu)PRO") 



Then for the Gilbert damping, since f = + F^ Aa changes as follows: 



Aa — > Aa 



1 - — In 



Mr 



(60) 



(61) 



d// (pR(y")j J PLiM)PR(M) 2 d? 

This result indicates that when writing data to an STT MRAM cell, the damping of the mag- 
netization dynamics is influenced by not only the spin pumping but also the bias voltage. 
However, the eff"ect of the bias voltage on Aa vanishes when both electrodes have the same 
electronic structure. 

In summary, we derived, at the microscopic level, the charge and spin currents passing 
through an MTJ in response to arbitrary motion of the magnetization in the free layer. The 
charge current consists of both Fermi surface and Fermi sea terms. The Fermi surface term is 
driven by the spin electric field and manifests as a dc current for steady precession of Ml(0 
in the direction of Mr, whereas the Fermi sea term is due to virtual transitions and essen- 
tially manifests as the ac current. With regard to spin transport, we focused particularly on 
the enhanced Gilbert damping (or the spin pumping effect) and thus obtained the microscopic 
expression for the enhanced Gilbert damping constant Aa = j^;^- Under a bias voltage, 
the DOSs of the two layers in the MTJ are shifted. Thus, the bias voltage changes the ef- 
fective exchange constant and the conductances, thus producing modulation of Aa. All the 
conductances consist of the tunneling amplitude Tlr and the local DOS on the surfaces of 
the layers; the real part of a retarded Green's function can be obtained from the imaginary 
part (namely, the local DOS) via the Kramers-Kronig relationship. In this formulation, the 
properties of the barrier layer material are considered in the local DOS, which can be easily 
obtained by first-principles calculations. 
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